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Abstract 

We present a method for the resolution of (oscillatory) nonlinear problems. It is based on the 
application of the Linear Delta Expansion to the Lindstedt-Poincare method. By applying it to 
the Duffing equation, we show that our method substantially improves the approximation given 
by the simple Lindstedt-Poincare method. 
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I. INTRODUCTION 



In this letter we present a method for the resolution of oscillatory nonlinear problems. 
This new method is obtained by applying the Linear Delta Expansion (LDE) l| to the 
well-known Lindstedt-Poincare (LP) method We explicitly solve the Duffing Equation 
(Anharmonic Oscillator Equation) and find that our approximation works much better over 
a wider range of parameters than does the simple LP method. We also compare our results 
to an approximation performed using the method of the perturbative S expansion ,3|, and 
find again that our approximation works with better accuracy and convergence for a wider 
range of parameters. 

In Section IH] we briefiy review the LP and LDE approximation methods, then in Sec- 
tion mil we show how by applying the LDE to the LP method, one can solve the Duffing 
Equation. Our results are presented in Section IIVI and we finally present our conclusions 
and directions for future work in Section IVl 



II. APPROXIMATION SCHEMES 



A. The Lindstedt-Poincare method 

In this section we introduce the Lindstedt-Poincare distorted time (LP) method We 
consider a nonlinear ODE of the form 

x(t) x{t) = e /(x(t)) , (2.1) 

which describes a conservative system, oscillating with an unknown period T. The nonlinear 
term e f{x{t)) is treated as a perturbation. Unfortunately, when the ordinary perturbation 
is applied to eq. 1)2.11) . by writing the solution as a series in e, the appearance of secular 
terms spoils the expansion and any predictive power is lost for sufficiently large time scales. 

In order to avoid the appearance of secular terms, we switch to a scaled time r = 27it/T = 
Q t, where T is the (unknown) period of the oscillations. The ODE now reads: 

n'^{T)+u;'xir)=ef{x{T)) . (2.2) 

We notice that the dependence upon e in this equation enters both in the solution x{t) 
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and in the frequency Q. By assuming e to be a small parameter we write 

and expand the r.h.s of eq. ()2.2|) as 
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By using these expansions inside eq. ()2.2|) we obtain a system of linear inhomogeneous 
differential equations, each corresponding to a different order in e. Let us consider the first 
few terms. To order we obtain the equation 
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describing a harmonic oscillator of frequency Q 
equation 
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where the r.h.s. is given by 
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uj. To order e we obtain the 
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We remark the oscillatory behavior of the driving term Si(r), because of its dependence 
upon the order-0 solution, xo(t). As a result si(r) will contain the fundamental frequency, 
corresponding to a period of 2tt in the scaled time, and multiples of this frequency, appearing 
through the term /(xo(t)). The presence of a driving term with the fundamental frequency 
leads to a resonant behavior of xi{t) and to the unfortunate occurrence of secular terms, 
which spoils our expansion. However, we can deal with this problem by fixing the coefficient 
ai to cancel the resonant term in the r.h.s. of eq. ()2.4j) . The iteration of this procedure to 
a given order n allows to determine the coefficients ao, . . . , a„ and therefore the frequency 

B. Linear delta expansion 

The linear delta expansion (LDE) is a powerful technique which has been originally 
introduced to deal with problems of strong coupling Quantum Field Theory, for which the 
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naive perturbative approach is not useful. Since then this method has been apphed to a 
wide class of problems jj, S, 0, S 1^ • In its original formulation a lagrangian density £, 
which is not exactly solvable, is interpolated with a solvable lagrangian Co{fJ,), depending 
upon one (or more) parameters fi: 

Cs = Coil^) + 6 (C-CoiiJ.)) . (2.6) 

For (5 = one obtains CoifJ^), whereas for S = 1 one recovers the full lagrangian Cs- The 
term 6 {C — Co) is treated as a perturbation and 6 is used to keep track of the perturbative 
order. Eventually 6 is set to be 1. 

We notice that the interpolation of the full lagrangian with the solvable one, Co{fi), brings 
an artificial dependence upon the arbitrary parameter fi. Such dependence, which would 
vanish if all perturbative orders were calculated, can be milden to a finite perturbative order, 
by requiring some physical observable O to be locally insensitive to /i, i.e: 

This condition is known as Principle of Minimal Sensitivity (PMS) and is normally seen to 
improve the convergence to the exact solution. 

III. ANHARMONIC OSCILLATOR 

Following the techniques described in the previous Section, we can now apply the LDE 
to the LP for the solution of the Duffing Equation (the Anharmonic oscillator equation): 

^it)+u^ x{t) = -fix\t) . (3.1) 

This equation describes a conservative system, where the total energy is given by 
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The period of the oscillation can be calculated in terms of an elliptic integral 



(3.2) 
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Texact = 2 / dx , (3.3) 

^2{E-V{x)) 



where A is the amplitude of the oscillations. 
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Following the procedure explained in the previous Section, we write Eq. fl3.ip as 



(3.4) 



where an arbitrary parameter A with dimension of frequency has been introduced. Clearly 
for 5 = 1, Eq. ()3.4j) reduces to Eq. (jH.lj) . We repeat the procedures previously explained and 
find a hierarchy of linear inhomogeneous differential equations to be solved sequentially. 
To zeroth order we obtain the equation 



+ K + A') xo{t) = , 



with solution 



(3.5) 



Xo(r) = A COST . 



(3.6) 



The zeroth order frequency is then given by 



(3.7) 



We proceed to compute the first order and find that 

«o^ + (^' + A^)xi(r) = S^{T) 



3AV' 



where 

5*1 (r) = A COST ai + y 

Now «! is fixed by eliminating the term proportional to cost: 

3AV ,2 

We obtain the solution 



cos3r . 



^i(^) = - oo/ ,9 , COST + ^ ^ cos3r , 



and the frequency 
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which is observed to be independent of A. 
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The next order gives: 
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where now 
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cos 5t . 



128 (a;2 + a^) 

As before a2 is fixed by ehminating the term proportional to cost: 

3^4 /i^ 



q;2 = 



128 (cu2 + A2) ■ 
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We obtain the solution 



(23>1V - 32A2) A3//(-3A2// + 4A2) 

3^2(7") = n . .ovn COST + 



+ 



1024(a;2 + p)2 



128(u;2 + p)2 



COS 3t 



1024(u;2 + ^2)5 



cos 5t 



and the frequency 



fl^ = ao + ai + a2 — ou'^ + 



(3.15) 



(3.16) 
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Note that at this order the frequency now depends on the arbitrary parameter A. However, 
due to the explicit dependence, by applying the PMS, we would obtain the same solution as 
in the simple LP method. In order to get a different solution, we must go to the next order 
in the expansion. 

Finally, following the same procedure, we obtain the expression for the third order: 
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By eliminating the term proportional to cost we determine to be 



and the solution 



/i 547 /i^ _ 1472 + 1024 A 
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"^^"^ = -32768 i}^T^^ '''' 

A^ 297 /i^ - 768 A^ fi + 512 A^ 
+ 16384 (A2 + cu2)3 ''"''^^ 

^5 (-3 /x + 4 A^) /i^ 1 

+ 2048 (A2 + o;2)3 + 32768 (A^ + u^y ' 

The frequency to order is now obtained to be 

^ = "° + + + = ^ - 128 (u^ + A^) + 512 (A^ • ^'-'"^ 

This time, the frequency depends upon the arbitrary parameter A in a nontrivial way and 
we can apply the PMS in order to fix the value of A. We do this by imposing that = 0, 
which leads to the following result: 

A = A^. (3.21) 

Notice that since A depends linearly upon A the formula for obtained in this case does 
not simply correspond to an expansion in A. As a matter of fact we find that the frequency 
corresponding to this value of A is 

^2 ^ 64 A^ /i^ + 192 A'fiu;' + 128 

96 /i + 128 ■ ^ ' ' 

Notice that the Duffing equation (|3.ip is left invariant under the simultaneous rescaling of 



the anharmonic coupling /i and of the amplitude, i.e. fx ^ n' and A A' = A ^ fx/ fx'. This 
invariance is manifest in the equation ()3.22|) . which is function of A"^ /j,, which is invariant 
under this rescaling. 



IV. RESULTS 



We now present the results following from this analysis. In Fig. 1 we compare the exact 
frequency, calculated with Eq. (j3.3p with the frequency obtained with our method (LPLDE), 
equation ()3.22|) . and with the LP method, equation ()3.20p taking A = 0, both to third order 
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in perturbation theory. We take uj = jj, = 1 and vary the amphtude of the oscillations. We 
observe that our method yields an excellent approximation to the exact result even for large 
amplitudes, where the simple LP approximation fails. 

In Fig. 2 we compare the period obtained with our method to the exact period of Eq. (j3.3p 
and to the one obtained with the formulae of P], which are obtained by applying the nonlin- 
ear delta expansion. Our method provides an excellent approximation to the exact period 
over a wide range of the parameter /i, which controls the nonlinearity. The plots are ob- 
tained assuming uo = 1 and the boundary conditions x(0) = 1 and x(0) = 0. The formulae 
of behave badly in the region /i < 0, which corresponds to a potential well of finite 
depth centered around x = 0, and yield a precision comparable to the one achieved with 
our method for yU > 0. Corresponding to the value = the oscillator is in a position of 
(unstable) equilibrium and the exact period diverges. Notice that for large values of /i all 
the methods seem to give a good approximation to the exact solution, including the LP 
method (to first order), which (to third order) was behaving poorly in the case previously 
studied. Unfortunately the equations of are not suitable to be analyzed as in Fig. 1, and 
thus we cannot fully test the efficiency of this method. 

In Fig. 3 we plot the relative error corresponding to the different approximations for 
/i > 0. Our method to third order in perturbation theory yields an error typically smaller 
than the errors of the other methods and with a magnitude of about 0.1 %. 

In Fig. 3 we plot the relative error corresponding to the different approximations for 
/X > 0. Our method to third order in perturbation theory yields an error typically smaller 
than the errors of the other methods and with a magnitude of about 0.1 %. 



V. CONCLUSIONS 



We have presented a method for the solution of (oscillatory) nonlinear problems. It is 
based on the application of the Linear Delta Expansion to the Lindstedt-Poincare method. 
We applied it to the Duffing Equation and find that the new method converges faster and 
with greater accuracy than the simple LP method. Also, by comparing it with methods 
based on the perturbative 5 expansion, we show that our solution not only converges faster 
and more accurately, but it also works for a much wider range of parameters. We are 
currently applying our method to a wider class of nonlinear problems and we are also 
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interested in considering its possible extension to quantum systems. 
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FIG. 1: Squared frequency of the anharmonic oscillator as a function of the amplitude (arbitrary 
units). Lo = II = 1. 
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FIG. 2: Period of the anharmonic oscillator. 
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FIG. 3: Error corresponding to the different approaches for the case studied in Fig. 3 
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